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We investigate the electronic and magnetic properties of single Fe, Co, and Ni atoms and clusters 
on monolayer graphene (MLG) on SiC(OOOl) by means of scanning tunneling microscopy (STM), X- 
ray absorption spectroscopy. X-ray magnetic circular dichroism (XMCD), and ab-initio calculations. 
STM reveals different adsorption sites for Ni and Co adatoms. XMCD proves Fe and Co adatoms 
to be paramagnetic and to exhibit an out-of-plane easy axis in agreement with theory. In contrast, 
we experimentally find a nonmagnetic ground state for Ni monomers while an increasing cluster size 
leads to sizeable magnetic moments. These observations are well reproduced by our calculations 
and reveal the importance of hybridization effects and intraatomic charge transfer for the properties 
of adatoms and clusters on MLG. 

PACS numbers: 73.20.Hb, 73.22.Pr, 71. 15. Mb, 78.70.Dm, 68.37.Ef 
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In the last decades, the miniaturization of devices has 
been a main goal in fundamental and applied research. 
Understanding the smallest magnetic units consisting of a 
few atoms or even a single atom is a crucial step towards 
this direction. While magnetic nanostructures on met- 
als [l|, |2| , semiconductors [3| and insulators [^ have al- 
ready been intensively studied, the discovery of graphene 
allowed to investigate a new substrate with intriguing 
electronic properties [5]. The electronic structure of 
monolayer graphene (MLG) exhibits a linear energy dis- 
persion with massless Dirac fermions which offers new 
opportunities for spintronic devices. Tuning the band 
structure of graphene and at the same time realizing the 
desired spin textures is a challenging task. For this reason 
intensive studies on the disorder in graphene introduced 
by vacancies [6[, adatoms [7|, and various molecules |8| 
have been performed in the past. On the other hand the 
nature of the interaction between graphene and individ- 
ual transition metal (TM) impurities and how the impuri- 
ties' properties can be modified has been debated mainly 
in theory (cf. Refs. |9l4l l\ and references therein) and 
only a few experiments have been reported so far [71, ll2| . 
Especially magnetic properties of single magnetic adsor- 
bates on graphene have experimentally not been investi- 
gated so far. 

In this article, we study individual 3d TM atoms on 
graphene prepared by thermal decomposition of the SiC 
surface in Ar-atmosphere [l3|. This preparation method 
results in high quality monolayer graphene, see the sup- 
plementary material [14] for further details. By means of 
X-ray magnetic circular dichroism (XMCD) we confirm 
the presence of sizeable spin and orbital moments and 



high-spin ground states of Fe and Co adatoms. Instead, 
Ni adatoms are found to be nonmagnetic but show un- 
expected X-ray absorption (XAS) spectra. While a van- 
ishing XMCD- signal suggests a 3d^^4:S^ configuration, in 



line with Ref. [l2|, XAS shows a multi-peak structure. 



Thus, XAS indicates a no n- fully occupied d-shell of the 
Ni adatoms. We will explain these results based on ab- 
initio calculations using density functional theory (DFT). 

The experiments were performed in two separate vac- 
uum chambers. The scanning tunneling microscopy 
(STM) measurements were performed inside ultrahigh 
vacuum using a low-temperature STM operating at 5 K. 
The XMCD measurements were performed at the ID08 
beamline of the European Synchrotron Radiation Facil- 
ity. Transition metals were deposited in minute amounts, 
i.e. 0.005 — 0.014 monolayer equivalent (MLE), at a sub- 
strate temperature of 10 K (in both systems) in order to 
obtain well-isolated monomers. The oxide free surface 
during the XMCD study was verified by the absence of 
a peak at the oxygen K edge in X-ray absorption spec- 
tra recorded before and after metal evaporation. XAS at 
the Ls^2 edges was performed in the total electron yield 
mode using circularly polarized light at 10 K. Fields of 
up to 5 T were applied to align the magnetization at 0° 
(normal) and 70° (grazing incident angle) with respect 
to the surface normal. All spectra have been background 
corrected. 

To gain a detailed understanding of the XAS of Ni 
adatoms we performed atomistic simulations of the sys- 
tem using first-principles DFT. Two different packages 
were used, namely the Vienna ab-initio simulations pack- 
age (VASP) [15| with a projector augmented plane wave 




FIG. 1: (color online) (a) and (b): STM topographies of pris- 
tine graphene. Tunneling parameters are U = 1 (—0.1) V 
and / = 0.5 (0.5) nA. (c): Ni on graphene. The inset 
reveals hollow site adsorption. Tunneling parameters are 
U = 0.4 (-0.1) V and / = 0.1 (0.15) nA. (d): Co on graphene. 
The inset exhibits top site adsorption. Tunneling parameters 
are U = 0.4 (-0.5) V and / = 0.025 (0.3) nA. 



(PAW) based code [Tel, [I^, and the Wien2k pack- 
age [l8[ with a full-potential linear augmented plane wave 
(LAPW) [l9[ basis sets. We employed the generalized 
gradient approximation (GGA) |20, _21] to the exchange- 
correlation potential and converged the results carefully 
with respect to the number of k-points. The system was 
modeled by a 3 x 3 graphene supercell of 18 carbon atoms 
with Ni monomers or Ni dimers placed on top. A vac- 
uum of at least 8 A (Wien2k) and 22 A (VASP) was in- 
cluded. X-ray absorption spectra were calculated within 
the dipole approximation using Fermi's golden rule. We 
neglected core- hole effects, but we found in test calcula- 
tions that these do not significantly change the calculated 
x-ray spectra. 

In order to identify the adsorption geometry we per- 
formed high-resolution STM measurements, see Fig. [H 
An overview of MLG prepared on the silicon-terminated 
SiC(OOOl) surface before adatom deposition is presented 
in Fig. [iji. A high-resolution image of the surface is 
shown in Fig. [TJd exhibiting honeycomb as well as buffer 
layer corrugation [22|. We note that the surface is 
fully covered by MLG with a very low defect density 
which both are required conditions for spatially averag- 
ing XMCD measurements in order to avoid buffer/double 
layer graphene or defect contributions [1^ . Figures [Th 
and d show constant-current maps of MLG covered by 
single Ni and Co adatoms, respectively. The insets of 
subfigures [It and d present the adsorption sites for Ni 



and Co. In both insets the bright lines are guidelines 
for the eyes indicating rows of hollow sites. A close in- 
spection shows that in case of Ni the crossing perfectly 
fits the center of the adatom (hollow site - in line with 
Ref. [12|) while the center of the Co adatom is found 
to fit a triangle spanned by surrounding crossings (top 
site). These results perfectly match the predicted ad- 
sorption sites found by calculations using a general gra- 
dient approximation functional with an on-site Coulomb 
potential of /7 = 4 eV [sl, [lH . 

In order to exclude aggregation processes within the 
XMCD measurements we evaluate the diffusion barriers 
of the adsorbed Ni, Co, and Fe monomers. From cal- 
culated binding energies, we can already extract an ex- 
pected trend for the diffusion energy among the different 
TM adatoms on MLG. In a simple approximation the 
diffusion barrier is given by the difference of the bind- 
ing energies of different adsorption sites [9|. Using this 
assumption the diffusion barrier is predicted to be the 
largest for Ni adatoms. It is especially larger than for Co 
adatoms which itself is larger than the diffusion barrier 
of Fe adatoms. Indeed, we revealed by means of the STM 
experiments Ni adatoms to be most stable in agreement 
with theory. In case of Ni we find an effective diffusion 
temperature of Td ^ 97 K. Co adatoms are found to be 
less stable (Td ^ 50 K) than Ni adatoms. In case of Fe it 
was not possible to determine the adsorption site as well 
as the effective diffusion temperature Td because they 
were manipulated by the STM tip. 

For exploring the magnetic properties of the adatoms 
on MLG XAS and XMCD were employed. Fig. [2K-c 
shows the XAS and XMCD spectra of individual Fe, 
Co, and Ni impurities on graphene. For all impuri- 
ties the XAS spectra reveal a multi-peak structure at 
the Ls edge. The observed multi-peaks are contrary to 
XAS measurements performed on bulk crystals |24J . |25| 
and on monomers adsorbed on metallic surfaces (see 
Refs. uJ, |26|, 1271 and references therein) which generally 
hints toward a weaker hybridization with the substrate in 
case of MLG as for the aforementioned metal substrates. 
We note that to our best knowledge XAS /XMCD exper- 
iments on TM monomers adsorbed on graphite have not 
been performed yet. The XMCD signals (see Figs. [2^ 
and b) demonstrate the existence of non-zero magnetic 
moments for Fe and Co while the spin and orbital mo- 
ments for Ni are zero, cf. Fig. [2t. In the following we 
will first discuss the results for Ni, compare them with 
ab-initio calculations, and then turn toward Fe and Co. 

For a low Ni coverage, most Ni adsor bates are sup- 
posed to be monomers. Naturally, the electronic configu- 
ration of free Ni atoms is Sd^As^. However, on graphene 
a promotion mechanism of two electrons from the 4s to- 
wards the 3(i-shell fills up the Ni d-shell [12|. Therefore, 
we expect Ni monomers adsorbed on MLG/SiC(0001) to 
show a Sd^^As^ configuration in line with previous re- 
sults [l2|. For this ground state, on the one hand, no 
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FIG. 2: (color online) XAS and XMCD measurements for (a) Fe (0.005 MLE), (b) Co (0.005 MLE), (c) Ni adatoms (0.007 
MLE), and (d) higher Ni coverage (0.014 MLE) at T = 10 K, B = 5 T. Upper and lower panels display L3,2 edges of XAS and 
XMCD spectra achieved with parallel (/x+) and antiparallel (/i-) alignment between the helicity of the incident beam and the 
magnetic field B, obtained at normal (0°) and grazing (70°) incident angle. The spectra have been normalized with respect to 
the incident beam intensity and the L3 preedge intensity. The insets in subfigures (a), (b), and (d) show the XMCD divided 
by the averaged XAS intensity for 0° and 70° , respectively. 



XMCD signal should be observed (see Fig. [2t). On the 
other hand, the XAS spectra are expected to show only 
step-like features because 2p -^ 3d transitions are for- 
bidden and 2p -^ 4s transitions would result in step-like 
features ^]. This expectation relies on the common as- 
sumption of a broad 4s band. In contrast, we observe a 
fine multi-peak structure for Ni adatoms, which usually 
is observed for non- fully occupied d-shells of Ni atoms, 
as in case of measurements performed for Ni in a 3d^ 
configuration in the gas phase [29|. 

To understand the multi-peak structure we simulated 
local density of states (LDOS) and XAS for Ni adatoms. 
The Ni atoms were placed above the hollow site where 
the relaxed height is 1.55 A. The theoretical monomer 
spectra (see Fig. ^) exhibit a multi-peak structure at 
the I/3 edge with two pronounced peaks, one at 0.7 eV 
(A') and a second at 1.6 eV (B'). We attribute these 
peaks to the two peaks in the experimental spectra at 
854.8 eV (A) and 856.1 eV (B), shown in Fig. Et- While 
the experimental peaks are separated by 1.3 eV, we find 
0.9 eV in our simulations. According to Fig.[3t^ the XAS 
(black) mainly follows the (i-LDOS (red). Although the 
empty s orbitals provide the largest peak at 0.7 eV, they 
show no significant contribution to the XAS because the 
transition matrix elements for 2p -^ 4s are much smaller 
than for 2p -^ 3d excitations. This is due to the nodal 
structure of the 4s wave function and the strong local- 
ization of the Ni d bands [30|. Hence, the peaks in the 
monomer case originate almost entirely from unoccupied 
Ni d orbitals. One possible scenario to explain the origin 
of the available unoccupied d-LDOS is that the electronic 
configuration of Ni is close to d^ or d^ and exhibits a 



magnetic moment. However, this is in contradiction to 
the experimentally determined absence of an XMCD sig- 



nal and also in disagreement with previous studies |12 |. 
Our DFT simulations confirm that the Ni monomer ex- 
hibits no magnetic moment. Hence, Ni is close to the 
d^^ configuration and therefore exhibits a much higher d 
occupation than bulk Ni. This causes the absence of a 
local magnetic moment in the Ni adatom: analogous to 
the case of an Anderson impurity model, there can be no 
magnetic moments formed close to full occupation. How- 
ever, a small amount of hole states is still present in the 
Ni d-shell which originates from a hybridization of the Ni 
d-orbitals with the Pz states of graphene. The hybridiza- 
tion is anisotropic which leads to the existence of multiple 
peaks in the LDOS and XAS. Our calculations find that 
the A' peak in Fig. [3^ is caused by the hybridization of 
Ni dz^ with Ni s orbitals which marginally contribute to 
the XAS. In contrast, the B' peak is dominated by a hy- 
bridization of d^y and d^2_y2 states with the graphene 
Pz states. This particular orbital character indeed traces 
back to the symmetry of the states in the conduction 
band near the van Hove singularity in graphene, which 
couple exclusively to this kind of d orbitals (cf. Ref. [23|). 
Increasing the coverage of Ni on MLG causes aggre- 
gation and therefore leads to clusters on the surface. 
Again, we find a strong multi-peak feature in the XAS 
spectra, but most importantly, we now observe a sig- 
nificant XMCD signal at the L3 edge, see the inset in 
Fig. EH. To explain these findings XAS and LDOS of 
Ni dimers (with atoms placed on hollow positions; re- 
laxed height: 1.60 A) were computed because more ad- 
sorbates containing at least two Ni atoms will be formed 



(cf. Fig. [Sb). The experimental XAS spectra show two 
main peaks at 854.4 eV (C) and 855.6 eV (D) (Fig.EJi). 
In agreement we find two main peaks at 0.2 eV (C) and 
1.4 eV (D') in our simulated spectra. Hence, the rel- 
ative intensities between (C) and (D) and the peak-to- 
peak distance are well-reproduced. Again, XAS mainly 
follows the Ni d-LDOS because the 2p -^ As transition 
matrix element is small. The C peak, corresponding to 
A' peak in the monomer case, is shifted downwards by 
^0.5 eV. This downshift can be attributed to changes 
of the (i-LDOS upon cluster formation. A comparison of 
the monomer and dimer LDOS shows that the unoccu- 
pied dimer s states are shifted towards the Fermi level 
while more d states above the Fermi level occur. Hence, 
s orbitals of both Ni atoms hybridize which partially pop- 
ulates s bonding states and depopulates the d orbitals. 
This effect amplifies for larger clusters. Therefore, size- 
able magnetic moments can be obtained for large clusters 
while Ni monomers are non-magnetic since they are close 
to a d^^ configuration. The Ni dimers turned out to be 
non-magnetic as well - from additional first-principles 
calculations we found that Ni clusters with more than 
four atoms are necessary for a finite magnetic moment. 
The magnetic properties of Ni clusters, crucially depend- 
ing on the number of atoms n, can thus be tuned be- 
tween non-magnetic {nm < 4) and magnetic {nm > 4) 
behaviour whereby magnetic moments of about 0.85 /iB 
per Ni atom can be achieved in the limit of infinitely large 
clusters, i.e. one monolayer of Ni adsorbed on graphene. 
In case of Fe we observe a pronounced multi-peak 
structure in the XAS spectra. The origin of this structure 
remains unclear and it might be traced back to manybody 
effects, crystal field, and hybrdization effects. The pre- 
dicted electronic configuration of Fe monomers on MLG 
is 36^^45^. The experimentally determined branching ra- 
tio (BR) [3l| has a value of 0.74 and, thus, indicates 
a high-spin ground state for the Fe monomers in agree- 
ment with theoretical predictions [QJ. M oreover, Fe has 
the smallest value for the ratio (R) [32| among all stud- 
ied TM, i.e. i?Fe = 0.05 ± 0.03, Rco = 0.35 ± 0.02, 
^Ni-ciusters = 0.42 ± 0.02. Hcucc, wc coucludc that Fe 
atoms exhibit a small orbital moment, in agreement with 
Ref. [33|. In addition, we determine the easy axis for 
the magnetic moments from the ratio of the XMCD and 
the XAS signal strength [27]. For Fe we find this ratio 
for normal incidence being 20 % larger than for grazing 
incidence indicating an out-of-plane easy axis (inset in 

Fig.EK). 

Concerning the XAS measurements for Co the multi- 
peak structure is less pronounced. We relate this to a 
stronger hybridization with MLG as compared to Fe. 
This observation is particularly in line with the pre- 
dicted binding energies of Co and Fe (0.62 eV and 
0.27 ± 0.02 eV, respectively) [R]. The branching ratio 
of ~ 0.74, extracted from the XAS spectra, indicates 
a high-spin ground state for Co as well and is in line 
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FIG. 3: (color online) Theoretical XAS spectra as well as the 
s- and d-LDOS of (a) Ni monomers and (b) Ni dimers on 
graphene. All spectra and LDOS are spin averaged. 



with a 36^^45^ electronic configuration of the Co adatoms. 
The analysis of the relative XMCD intensity of normal 
and grazing incident reveals an out-of-plane easy axis in 
agreement with Ref. [33|. 

In summary, we find top site adsorption for Co 
monomers while the adsorption site of Fe adatoms could 
not be determined due to tip-induced atomic manipula- 
tion during our STM measurements. For Fe and Co we 
observe a multi-peak structure in XAS and a non-zero 
XMCD signal. Both species are paramagnetic and ex- 
hibit an out-of-plane easy axis while the branching ratio 
indicates high-spin ground states. For single Ni atoms we 
observe a multi-peak structure in XAS as well while the 
XMCD shows a zero magnetic moment. Our ab-initio 
calculations show that this multi-peak structure is due 
to a small amount of holes in the almost entirely full 
Ni 3(i-shell. These holes are due to hybridization of Ni 
and graphene. Further Ni deposition results in clustering 
and the rise of non-zero magnetic moments. The main 
features of the experimental XAS spectra are well repro- 
duced by simulating Ni dimers. Moreover, first-principles 
calculations reveal a transition from non-magnetic clus- 
ters toward magnetic clusters when the Ni cluster size 
exceeds four atoms, resulting in sizeable magnetic mo- 
ments. 
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The particular samples used in this work were prepared by graphitization of the Si-face, 
the so-called SiC(OOOl) surface, in Ar-atmosphere fll]. As a consequence of this method 
almost defect free carbon layers can be achieved. Figure ISTk presents a large scale AFM 
image of such a sample. The average terrace width of our samples is found to be about 
300 nm. The observable step edges are subsurface steps of the underlying SiC substrate. 
Figures ISTb and c show a STM constant current map of a step edge and the respective line 
profile. The measured apparent height of 2.5 A matches the spacing of one SiC bilayer. 
The carbon layer itself is an uninterrupted and closed film continuing across the SiC step 
edge j2| . Neither regions of double layer graphene nor regions of an uncovered buffer layer 
can be found which indicates that the desired coverage of one monolayer of graphene was 
perfectly achieved. 
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FIG. SI: (color online) (a): Large scale topography of MLG on SiC(OOOl) measured by means 
of atomic force microscopy in tapping mode, (b) STM topography of a subsurface SiC step edge 
indicating the absence of carbon buffer layer or bilayer graphene areas. Tunneling parameters are 
(7 = — 1 V and / = 0.5 nA (c) Line profile of the SiC step edge as indicated in (b). 



After the low temperature deposition of Ni or Co atoms onto MLG/SiC(0001) we did not 
find any attractive or repulsive interaction between the adatoms and the SiC step edges by 
means of scanning tunneling microscopy. The adatoms are statistically distributed across 
the entire surface. Neither the local density of states (LDOS) of adatoms nor the LDOS 



of the graphene layer is varying near these subsurface step edges which has been proven 
by scanning tunnehng spectroscopy. Besides the experiments performed directly after low 
temperature deposition, we heated the samples to variable temperatures above and below 
the diffusion barrier. We note that single atoms do not interact with the step edges, i.e. the 
step edges do not change the surface potential of graphene and do not modify the properties 
of the adatoms. Hence, we exclude step edges to cause multi-peak XAS spectra. 

Lauffer and co-workers showed that graphene is bent in the vicinity of an underlying step 
edge J2|. This causes an indirect influence on the adatoms, e.g. the magnetic anisotropy di- 
rection of atoms adsorbed in the bent region of graphene is tilted compared to the anisotropy 
direction of monomers adsorbed on a flat terrace. If we assume an area of 2 nm (eight unit 
cells of graphene) in which graphene is not flat, we can estimate the number of adatoms 
adsorbed in this region from the average terrace width. We expect 0.6 % of all adatoms 
to adsorb on bent graphene, resulting in a negligible influence on the measured XAS and 
XMCD spectra. 
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